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ABSTRACT 


This study was undertaken to determine if the US Army's M119 
Howitzer could be safely towed in the “firing” position using the 
existing A-Frame brace to support the gun tube during movement. A 
computer model of the howitzer was developed that would predict 
motions of the howitzer and consequently, loads on _ critical 
components. The results obtained from this computer model were 
validated by comparing them to actual data obtained from = an 
instrumented test done on the howitzer at Aberdeen’ Proving 
Grounds, Md. The validated model was used to predict the "worst 
case" loads on the  howitzer's Roseme. and the T-Bar _ during 
movement in the "firing" position. The maximum predicted stresses 
in the A-Frame and the T-Bar from the “worst case" loads were 
compared to the yield stress. These comparisons showed that the A- 
Frame and T-Bar were strong enough to withstand the “worst case" 
predicted loads that would be caused by towing the M119 Howitzer 


in the "firing" position. 
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I. INTRODUCTION 


A. GENERAL 

The purpose of this study was to determine the feasibility of 
transporting the M119 howitzer in the firing position over terrain 
and at speeds used during normal tactical operations. The’ end 
objective was to determine if costly design changes would be 
required to safely permit towing the M119 MHowitzer in the firing 
position. The approach used was to model the vertical plane 
dynamics of the towing vehicle (HUMMV) and the howitzer, model 
representative terrain inputs to the wheels in the form of potholes, 
and to perform experiments to validate the model. It is proven that 
the simulation model is reasonable and it follows that stresses in the 


critical components are not expected to warrant redesign. 


B. BACKGROUND 

The M119 Howitzer is a British designed howitzer that is 
currently being built at Rock Island Arsenal (RIA), [Illinois and 
fielded to US Army Field Artillery units. The M119 Howitzer 
Operational manuals specify towing procedures that require the gun 
tube to be placed in the stowed position when moving the howitzer 
during tactical operations. These towing procedures are’ different 
and more difficult to implement than those used for all other towed 
howitzers in the US Army. All other US Army towed howitzers are 


towed during tactical operations with the gun tube in the firing 


position, facilitating rapid emplacements and displacements that are 
required for successful tactical operations. 

The M119 Howitzer has an existing A-Frame brace which can 
Support the gun tube when the howitzer is towed tin _ the firing 
position. The primary concern of this study was whether or not the 
A-Frame and the T-Bar structure on the howitzer's trail section, to 
which the bottom of the A-Frame brace connects, will be strong 
enough to support the loads generated when towing the howitzer 
over rough terrain. This study also addresses concerns about trail 
loads on the howitzer. This study does not evaluate rollover stability 
for the M119 Howitzer when towing in the firing position. 

The terrain that the HUMMV/Howitzer system can traverse is 
Clearly limited by the terrain the HUMMV can negotiate, therefore, 
the Mil9 Howitzer and HUMMV have been evaluated as A system. A 
dynamic simulation model of this system was seen to be a _ necessary 
tool to aid in making design decisions from a rational viewpoint and 


is the basis of this study. 


C. METHOD OF ANALYSIS 

To analyze the M119 Howitzer a computer model of _ the 
HUMMV/Howitzer system was developed to predict the response of 
the system to roadway inputs. This computer model treats the major 
mass elements of the HUMMV/Howitzer system as lumped masses to 
develop vertical plane heave and pitch equations of motion for the 
system. The development of this computer model is contained in 


Chapter II. To drive the computer simulation model, terrain models 


were developed for both a "worst case" trench type pothole and a 
sinusoidal "washboard" roadway. Details on the development of 
these roadway inputs are contained in Chapter III. Chapter IV 
contains the results of driving the computer model over the roadway 
inputs. Chapter V explains the structural analysis that was used to 
determine the maximum stress in the A-Frame and the T-Bar. 
Chapter V also contains the procedures and results of calibrating an 
A-Frame that was used as a load cell for an instrumented dynamic 
test. Chapter VI explains the conduct and the results of the 
instrumented test on the M119 Howitzer/HUMMV system that was 
conducted at Aberdeen Proving Grounds, Maryland in April 1993. 
Also, in Chapter VI the results of the model are compared to _ the 
experimental results from the instrumented test. This comparison 
showed reasonable correlation between the model and experimental 
results. Conclusions and recommendations based on this study are 
contained Chapter VII. The final recommendation is that towing in 


the firing position should be allowed. 


Il. TUMMVIHOWITZER SYSTEM COMPUTER MODEL 


The computer model developed is based only on the heave and 
pitch equations of motion (EOM) for the HUMMV/Howitzer system. 
The reason for this simplification is that the lateral, yaw, and roll 
motions are much less significant for determining the forces in the 
A-Frame for the type terrain inputs considered by this study. While 
horizontal plane motions are key elements for the analysis of turns 
and rolling during transit, the concerns of this study indicate that a 
vertical plane model is adequate. The major mass groups of the 
HUMMV/Howitzer were’ treated as lumped ~omasses_ in _ the 
development of the equations of motion. The equations of motion 


were developed using Lagrange techniques. 


A. HUMMV/IHOWITZER SYSTEM PARAMETERS 

The primary sources of information for system parameters for 
the HUMMV and _ howitzer were the US Army _ Tank-Automotive 
Command Research, Development & Engineering Center Technical 
Report No. 13337, [Ref. 1], for HUMMV data and US Army Combat 
Systems Test Activity, APG, Report No. P-83179, [Ref. 2], for the 
M119 Howitzer data. Additionally, many system parameters .were 
determined from direct measurement on the HUMMV and howitzer 
or from design diagrams, or obtained from the study sponsor at Rock 
Island Arsenal (RIA). For readers of this study that are unfamiliar 


with the components of the M119 Howitzer, Appendix G_ contains 


Several diagrams that depict the howitzer components’ that are 
important for this study. 

The system parameters used for the HUMMV are based on using 
the Model No. M1069 HUMMV as the prime mover for the M119 
Howitzer. The cargo loading of the HUMMV significantly effects its 
center of gravity and consequently effects all of the HUMMV's 
System parameters which are dependent on its center of gravity. 
Modeling was therefore limited to using the same _ cargo loading 
configuration that was used in’ [Ref. 1] and other loading effects were 
ignored. These HUMMYV simplifications are justified because for this 
computer model the HUMMV serves primarily as the driving force at 
the howitzer hitch causing pitching motions of the entire howitzer 
and we are not concerned with the loading of HUMMV components. 
The major components of the HUMMV/Howitzer system computer 
model are identified in Figure 2-1 below. 

In the model, the tires act as linear springs which provide force 
inputs to the HUMMV and howitzer wheels, based on the terrain over 
which the model was "driven". The suspension systems for both the 
HUMMV and the howitzer are modeled as_ springs and dampers 
connected in parallel. The A-Frame and T-Bar are modeled as a set 
of stiff springs connected in _ series. The following paragraphs 
provide more detail on the modeling of the individual components 
and the determination of appropriate parameter values. Appendix A 
provides a listing of the numerical values, units, and the symbols 
used for each system parameter in the computer model. The 


parameter values for those components with identical components on 
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Figure 2-1 


the right and left sides of the system (wheels and _— suspension 
components) are multiplied by two. The values in Appendix A_ for 
these components reflect the effects of both sides. 

1. HUMMV and Howitzer Tires 

The HUMMV has the same tires on the front and_ “rear, 
however, the front and rear tire pressure is different (20 psi vs 30 
psi) resulting in the rear tire spring being much _ stiffer. The 
stiffness used for the HUMMV tires was obtained from normal force 
vs tire deflection data [Ref. 1]. The stiffness of the howitzer tire used 
in the model was estimated because no measured data on the 
howitzer tire stiffness was available. 

The tires are modeled as linear springs which transmit a 
vertical force to the respective wheel masses proportional to the 
amount that the springs are stretched or compressed. The amount of 
force is determined by the relative dynamic displacement between 
the terrain input for the wheel and the center of mass of the wheel. 
If the tire leaves contact with the ground there will be no input from 
the terrain until it again returns to contact. To account for this loss 
of contact situation in the computer model, the tire spring force is set 
to the static weight that the tire supports when the tensile (de- 
compressive) spring force exceeds the static weight on the tire. The 
tire springs are turned back on in the computer model when the 
conditions are met for road surface contact to be re-established. 

Pap HUMMV Suspension System 
The HUMMV suspension. system consists of a _ shock 


absorber and coil spring for each wheel. The shock absorber passes 


through the middle of the coil spring. . The force in the spring is 
linearly proportional to the relative displacement between the center 
of mass of the wheel and the attachment point on the HUMMV body. 
The values for the spring constants were obtained from Table 5-13 
of [Ref. 1]. Precise modeling of the shock absorbers is_ difficult 
because they have a mid-stroke region and a hydraulic bump stop 
region. In the mid-stroke region the shock absorber force is linearly 
dependent on the relative velocity of the center of mass of the wheel 
and the attachment point on the HUMMV. However, the linear 
damping constant changes depending on whether the force in _ the 
Shock absorber is _ tensile or compressive. When the shock absorber 
is compressed or stretched into the hydraulic bump stop region the 
damping of the shock absorber greatly increases. This sudden 
increase in damping causes problems in the integration routine that 
was used for the computer model. To overcome these problems and 
still account for the energy loss due to the increased damping in the 
hydraulic bump stop region, the mid-stroke linear damping constant 
Obtained from [Ref. 1] was increased by 20% for both the front and 
rear shock HUMMV absorbers. Figures 2-2 and 2-3 shows the force 
versus velocity plots that were obtained from [Ref. 1] for the front 
and rear shock absorbers, respectively. 
3. Howitzer Suspension System 

Modeling the howitzer suspension system iS more 
complicated than the HUMMV because it consists of a torsional spring 
and a shock absorber that is not vertically oriented. The torsional 


Spring is modeled as a linear vertical spring that is connected to the 
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Figure 2-2 Force vs Velocity Plot for the Front HUMMV Shock 
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Figure 2-3 Force vs Velocity Plot for the Rear HUMMV Shock 
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center of mass of the howitzer wheel and attached to the howitzer 
trail directly above the wheel center of mass. The effective spring 
stiffness constant was calculated based on the measured data of the 
torsional spring stiffness. Figures 2-4 through 2-7 show the angular 
deflection versus vertical force plots that were contained in [Ref. 2]. 
Using the average torsional spring stiffness from this data and_ the 
length of the road arm, an angular rotation was converted to an 
equivalent vertical displacement at the end of the road arm. The 
spring constant for the model was then calculated by determining 
the constant value that must be multiplied times this vertical 
displacement at the end of the road arm, to produce a moment that is 
equivalent to that caused by the torsional spring for a_ specified 
angular rotation. 

A similar procedure was used to model the howitzer's 
shock absorber as a vertical shock absorber with a linear damping 
constant that is velocity dependent. The attachment points for the 
ends of the shock absorbers in the model are the same as_ the 
modeled end points of the howitzer suspension spring. The actual 
non-vertical shock absorber force was transformed to a vertical force 
in the model by determining the amount of vertical force applied at 
the end of the road arm needed to produce an equivalent moment 
about the road arm pivot point. Like the HUMMV shock absorbers, 
the howitzer shock absorbers are non-linear and have different 
damping values for tension and compression that are accounted for 
in the computer model. This non-linear damping in_ the shock 


absorbers is a significant aspect of the computer model. The _ values 


1 1 


used for the shock absorber damping constants were determined 
from force versus velocity test data provided by RIA. Figure 2-8 


shows the force versus velocity plot for the howitzer shock absorber. 
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Figure 2-4 Torsional Spring Stiffness for Serial No. 22 (Upward 
Motion) 
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Figure 2-8 Force versus Velocity Plot for Howitzer Shock 
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4. A-Frame/T-Bar Stiffness 
The A-Frame and T-Bar are modeled as_ stiff springs 
connected in series to determine a combined linear spring stiffness. 
The stiffness of the A-Frame was determined during the A-Frame 
calibration test (explained in Chapter V). The stiffness of the T-Bar 
was determined by using the principle of superposition of forces to 
determine the effective bending stiffness for the non-symmetric 
bending of the T-Bar that is caused by the A-Frame. This non- 
Symmetric bending of the T-Bar results because the angle at which 
the A-Frame transmits forces and moments to the T-Bar is _ not 
perpendicular to an axis of symmetry of the T-Bar. 
5. Equilibrator Springs 
The equilibrator springs provide only small dynamic forces 
because their combined stiffness is nearly 1000 times ess than the 
effective stiffness of the A-Frame/T-Bar assembly. Therefore, 
equilibrator springs are neglected in the computer model. 
6. Elevating Clutch 
The effects of the elevating clutch of the howitzer were 
neglected in the final computer model. It is believed backlash in the 
elevating mechanism gears will prevent the clutch from _ being 
activated while towing with the A-Frame attached because the 
relative rotation of the gun tube is very small. However, because of 
some uncertainty about the clutch several computer runs were made 
with clutch modeled as a linear torsional damper. The value of the 
torsional damper used in the model was estimated because no test 


data on the clutch was available. These model runs with the clutch 


Ni 


included, produced only slightly different results from the runs 
without the clutch included. Therefore, because of the previously 
mentioned belief that the clutch would not be activated, it was 
neglected in the final model. 
fhe HUMMV Component Mass Properties 

The HUMMV is divided into three lumped masses, _ the 
HUMMV body mass and two wheel masses. The front and _ back 
HUMMV wheels are treated as having identical mass values in the 
model. Each wheel weight/mass includes the weight of the tire, 
wheel, shocks, suspension springs, and suspension support arms. The 
HUMMV body weight/mass varies dependent upon it's cargo load. 
The location of the center of gravity also can vary and is determined 
by the positioning of the cargo. Therefore, the values used in _ the 
computer model for the HUMMV mass, the mass moment of inertia 
for pitch, and lengths that are based dependent on the center of 
gravity, are estimated values. These estimated values are based on 
the same loading situation as was used in [Ref. 1]. The pitch mass 
moment of inertia for the HUMMV is about the center of gravity of 
the sprung HUMMV mass. 

8. Howitzer Component Mass Properties 

The howitzer is divided into three separate lumped masses 
in the model; the howitzer wheel masses, the elevating mass, and 
the combined trail/saddle mass. The wheel weight/mass includes 
the tire, rim, and approximately 30% of the road arm that rests on 
the wheel assembly. The elevating mass includes the gun tube, the 


cradle, the recoil system and all other parts that move the when the 
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gun is elevated about the trunion. The mass moment of inertia of the 
gun tube for rotation about the trunion was estimated using the 
approximate masses of all the elevating components. Thies center of 
gravity location for the elevating mass was estimated by RIA. The 
combined trail/saddle mass includes all of the gun mass except the 
elevating mass and the wheel mass. The center of gravity location 
for the trail/saddle mass and the pitching mass moment of inertia for 
the trail/saddle mass was estimated from design diagrams. 
9, Length Properties 

All lengths used in the computer model are identified in 
Figure 2-9 below. These length properties were determined either 
by direct measurements on the howitzer or by scaling from design 


diagrams. 


B. EQUATIONS OF MOTION (EOM) 

In this section the heave and pitch equations of motion (EOM) 
for the model are developed from Lagrange's equations using the 
HUMMV/Howitzer system parameters that were previously 
explained. Lagrange's equations require the determination of the 
System's total kinetic and potential energy and the power dissipated 
by the dampers, in terms of a set of generalized coordinates. 

Throughout this development of the EOM, expressions are 
developed that are functions of generalized positions and _ velocities, 
angular displacements, angular velocities, and angular accelerations. 
To simplify the EOM with minimal loss of representational accuracy, 


these expressions are linearized using small angle approximations. 
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Figure 2-9 Length Properties of HUMMV/Howitzer System 
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The sine of a small angle is approximated as just the angle (sin0=6) 
and the cosine of a small angle is approximated as equal to one 
(cos0=1). These approximations are valid because the pitch or 
angular rotations are predicted to be less than five degrees. Order of 
magnitude considerations are used to further simplify the EOM. 
Because the angular displacements are small their angular velocities 
and angular accelerations are also small. Therefore, any terms that 
contain an angular displacement, angular’ velocity, or angular 
acceleration multiplied times another angular displacement, angular 
velocity, or angular acceleration are neglected because these’ terms 
will be very small. 


The general form of Lagrange's equations used is: 


d|{ OT OT 0U OF 
dt \ 0q; Oq; °q,; 94, 
where: 
Gi = generalized coordinates 
T = total kinetic energy of the system 
U = change in the potential energy of the system with respect to 


its potential energy in the static-equilibrium position 
F = the power dissipated by forces in dampers 
Q = the input force acting on the system. 


The development of the above components of Lagrange's equation 
and the assembly of these components into the system's heave and 
pitch EOM follows. Due to the length of this development the reader 
is provided only the final results from each step. The intermediate 


multiplication's, differentiation's, and linearization's are omitted. 


75) 


1. Generalized Coordinates (qj) 
The set of generalized coordinates, qj, that was used to 
_uniquely specify. the motion of the HUMMV/Howitzer system are 


depicted on Figure 2-10 and are defined follows: 


* Z,&Zo: The vertical displacements for the front and_ rear 
HUMMV wheel masses, respectively. 


* Z3: The vertical displacement of the center of gravity of the 
HUMMV body. 


* Z4: The vertical displacements of the center of gravity of the 
howitzer wheel mass. 


* 6,: The pitch angle of the HUMMV body about its center of 
gravity. 


* @5: The pitch angle of the howitzer trail/saddle mass about its 
center of gravity. 


* 63: The absolute rotation of the gun tube. 


As depicted on Figure 2-10, upward vertical motions and clockwise 
rotational motions are defined as positive. The equilibrium position 
for each of the generalized coordinates is the static equilibrium 
position for the respective components that are correlated to each of 
the generalized coordinates. In defining positions to key locations in 
the system unit vector notation of i, j, and k will be used 
corresponding to the x, y, and z axis directions of the Cartesian 
coordinate reference system shown on Figure 2-10. 

The kinetic and potential energy of the system and_ the 


damping and input forces are all computed in terms of these 
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2-10 Generalized Coordinates 


Figure 


System 


a, 


generalized coordinates and the road inputs U;,, U2, and U3. To carry 
out the required computations, position vectors to key locations on 
the HUMMV and howitzer are defined in terms of the generalized 
coordinates. The position vectors to these key locations are depicted 


on Figure 2-11 and defined as follows: 


* ry & rz = the front and rear attachment points, respectively, for 
the suspension springs and shock absorbers on the HUMMV. 


r3= the center of mass of the trail/saddle mass 


* rqg= the modeled attachment location of the suspension spring 
and shock absorber of the howitzer 


* rg= the location where the A-Frame attaches to the T-Bar 


* reo= the location of the center of mass of the gun tube/elevating 
mass 


* r7= the location where the A-Frame attaches to the gun 
tube/elevating mass 


These key locations that are defined above can be expressed as 
position vectors in terms of the above defined lengths (Figure 2-9), 
generalized coordinates (Figure 2-10), and the unit vectors 1, Jj, and k 


of the Cartesian coordinate reference system (Figure 2-10) as follows: 


=(Z, - L,sin(6,))k + (L,cos(6,))i 


+ L,sin(6,))k - (L, cos(6, ))i 


Z, 
Z, + 1, sinG ys L, ,sin(6,))k = ie cos(6,) + L, ,cos(8, ) i 


Z, + L,sin(0,) + L,,sin(@,))k - (L,cox@,) + L,,cos(0,))i 
=(Z, + L,sin(6,) + (L,,+L,,)sin(®,))k - (L,cos(6,) + (L,, +L,,)cos(@,))i 


( 
( 
( 
( 
( 


ry 
r, 
* r, 
ry, 
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Key Locations for the HUMMV/Howitzer System 
2D 


Figure 2-11 


L,cox8,) + L,,cos(6,)). 
+e. cos( Ott Gy) J 
L, cos(8,) + L,,cos(6,)). 
+L,.cos(@, + 0;) J 


*r =(Z, + L,sin(0,) + L,,sin(6,) + L,, sin(®, + 0,))kK - ( 


*r =(Z, + L,sin(0,)+ L,,sin(6,) + L, sin(0, + 0,))k — ( 


The velocities of these key locations can be determined by 
taking the time derivatives of the position vectors. However, only 
the velocities of locations where there are velocity dependent actions 
occurring in  Lagrange'’s equations need to be determined. The 
required velocities of the key locations on the HUMMV and_ howitzer 


are as follows: 


ll 


(Z, - L,cos(6,)6,)k - (L,sin(®,)6, )i 
=(Z, + L,cos(0,)6, )k + (L, sin(6,)6, 


( 
VA 4 L, cos(6, )0, iH L, ,cos(0,)6, )k + (ie sin(0, 0, + L, , sin(6,)6, ji 
(Z, + L,cos(6,)6, + L,,cos(0,)6,)k + (L,sin(6,)6, + L,,sin(6,)6, ji 
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In addition to the above velocities, the vertical velocity of the CG of 
the HUMMV, bee and each of the wheel masses, yA De. & Za, are 
needed in Lagrange's equation. 
2. Kinetic Energy (T) 
The velocities defined above are used to calculate the 
kinetic energy of the system. The kinetic energy is calculated by 


Summing the linear and rotational kinetic energy for each of the 
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system's masses. The total kinetic energy of the system is given by: 


¥ 1,=¥ (’my? + 4187) i= 16; j= 13 


where: 

Tj; = the kinetic energy of each mass in the system 

m;= the mass of each lumped mass 

vi = the linear velocity of the each mass 

I;= the heave mass moment of inertia of each rotating mass 
about its center of gravity 

Q0;= the angular velocity of each rotating mass 


Therefore, the total kinetic energy of the HUMMV/howitzer system 
was determined by summing the _ kinetic energy components as 


follows: 


T=4M,(Z,)’ + 4M,(Z,)’ + 4M,(Z,)’>+ $M,(Z,)? + 4M,(f3)” 
LM,(f)° + $1,(6,)° + 41,(6,Y + 41,(6,) 


In the above expression the square of the velocities is determined by 
Substituting in the appropriate velocity expression as defined above, 
and taking the dot product of the velocity with itself. This 
substitution and dot product operation is left for the reader. When 
taking this dot product all of the k dot i and i dot k cross _ product 
terms drop out because the 1 and k _ unit vectors are orthogonal. 
Having determined the total kinetic energy, the time derivative of 


the partial derivative of the total kinetic energy with respect to the 


; OT 
velocity of each of the generalized coordinates, FBG. can be 
determined. Carrying out this differentiation and linearizing yields 


the following: 


Za 


d OT 
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ars 36, = (M, ) ( *) 
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+ (1, + Moog); 

The total kinetic energy, T, must also be _ differentiated 
with respect to each of the generalized positional coordinates, 
however, for this system this differentiation is equal to zero for all 
coordinates. 

3; Potential Energy (U) 

The potential energy for this system results from _ the 
Strain energy that is developed in the suspension springs and the A- 
Frame spring. The gravitational potential energy that arises from 
changes in vertical positions is ignored because the resulting terms 
drop out of the dynamic equations of motion. The potential energy 


in any spring is determined as follows: 
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U, = 4K,(Ax,)’ 


i 


where: 
U; = the potential strain energy in each spring 
K; = the spring stiffness for spring 
Ax, = the change in length of the spring. 


Since locations to the attachment points of the ends of all the springs 
have been defined above in a common inertial coordinate’ reference 
system, taking the difference of the position vectors of the end points 
of each spring will yield the change in length of the spring. Thus for 
this system, the total potential energy is found by summing the 


potential energy in each of the springs as follows: 


Therefore, the partial derivative of the total potential energy with 


0U _ 
respect to each of the generalized coordinates, ——, after linearizing 


results in the following: 
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* = = -K,L,L,,cox6,)@, + K,L,L,,cos(6,)0; 


4, Generalized Power Lost Due to Damping Forces (F) 
The amount of generalized power that is dissipated by 


each of the dampers in the system given by the equation: 


F, = 4B,(Ax,)° 
where: 
F; = the generalized power dissipated in each damper 
B; = the damping constant for each damper 
Ax, = the difference in velocity of the attaching points of each 
damper. 


The total amount of generalized power dissipated in the system is 
found by summing the energy dissipated in each dampers. Sum ming 


the energy dissipated by the damping forces yields the following: 
a ‘ se \2 1 : SN? 1 , - \2 , - \2 
F =4B(r, = Zs) ? +B, (Fr, a Za) a +B,(i, - ZA) a 1B, (6, a 6, 


In Lagrange's equation this total dissipated energy, F, must’ be 
differentiated with respect to the velocities of each of the generalized 


OF 
coordinates, -——. 
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Carrying out this differentiation and again linearizing results in the 


following: 
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5. External Forces (Q;) 

The external forces acting on the system are the roadway 
inputs that are transmitted through the tire springs. For the 
computer model these input forces are treated as_ purely vertical 
forces. The forces result from the road inputs causing a change in 
length of the tire springs. These road inputs are depicted on Figure 


2-10 and defined as follows : 
* U; = the road input for the front HUMMV wheel 


* U>= the road input for the rear HUMMV wheel 
* U3= the road input for the howitzer wheel. 
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6. Assembling Lagrange's Equation 
Combining the parts of Lagrange's equations that were 
developed in the above paragraphs results in a system of equations 


of motion in the form: 


[M]{x} + [C(x)|{x} + [K]{x} = [f] 


where: 

[M] = the 7 x 7 system mass matrix 

[(C(x)] = the 7 x 7 system damping matrix that is a function of 
whether the shock absorbers are in tension or compression 

[K] = the 7 x 7 system stiffness matrix 

{x}, {x}, & {x} = 7 x 1 column vectors that contain accelerations, 
velocities and positions of all coordinates, respectively. 

[f] = the 7x1 forcing vector that results from the road inputs. 


The [M], [C(x)], and [K] matrices for this system are contained in 
Appendix B. This completes the development of the heave and 

pitch EQOMs. The resulting EOMs were numerically integrated to 
determine the time history for motions of the component masses of 


the system. 


C. HUMMV/IHOWITZER MODEL NUMERICAL INTEGRATION 
The following section describes the procedures that were used 
to numerically integrate the previously developed equations of 
motion. 
i Transformation to Ist Order Equations 
In order to numerically solve the equations of motion, the 
formulated system of second order ordinary differential equations 


must be transformed into a system of first order differential 


eZ 


equations. This equation transformation was performed using matrix 
algebra and results in the following expression for the system of first 


order differential equations: 


[ [0] [1] fx] fo ] 
[xdot] = - ee ices 
L-mky (=m c(or | Le” Lae] 
where: 
[xdot] = a 14 x 1 column vector of which the first seven elements 


are the velocities of all the generalized coordinates and the last 
seven elements are the accelerations of the generalized coordinates. 


[x ] 


lk | = a 14 x1 column vector of which the first seven elements 


are the positions in terms of the generalized coordinates and the last 
seven elements are the velocities of the generalized coordinates. 


{O] = a 7 x 7 zero matrix 
[1] = a 7 x 7 identity matrix 


[-M“'C(x)] & [-M"'K] = these are 7 x 7 matrices that result from 


multiplying the negative inverse of the system mass matrix, [M], 
times the damping matrix, [C], and the stiffness matrix, [K], 
respectively 


[oO | 


[Mf | = a 14 x 1 vector which the first seven elements are zero and 


the last seven elements are the matrix product of multiplying the inverse 
of the mass matrix, [MJ], times the forcing vector, [f]. 


For the integration the system was started from rest, thus, all 
seven of the system's displacement, {x}, and velocity, {x}, state variables 


had a value of zero for initial conditions. 
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| 
| 
2. Integration Difficulties ' 

The equations of motion were numerically integrated using the 
Matlab ODE 45 integration routine. This is a variable time step 4th order 
Runge-Kutta integration routine. Integrating these equations of motion | 
using ODE 45 was complicated by the non-linear damping of the _ shock 
absorbers and the required shutoff of the tire springs when they leave 
contact with the road surface. 

When the damping values of the shock absorbers would abruptly 
change from compression to tension and vice versa, there would be a 
corresponding abrupt change in the shock absorber force. Determining 
when the shock absorber damping value should change is built into the 
integration code by checking for each time step if each of the shock 
absorber forces are compressive or tensile. An "if" statement ensures that 
the correct damping value is used for the condition that exists. The abrupt 
changes in force in the shock absorbers slowed down the tolerance based 
ODE 45 routine. A_ similar situation occurred when the tire spring forces 
would abruptly change to the respective static weights when contact with 
the road was lost for a particular wheel. 

To increase the speed of the ODE 45 routine the integration was 
divided into one second time segments so the solution matrix that contains 
the state variables remains. relatively small. The displacements and 
velocity states at the end of the previous time segment are used as _ the 
initial conditions for the next one second time segment. 

The resulting time histories for the displacements and_ velocities 


allowed post-processing to be completed determine accelerations, forces, 
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and stress in key system components. The computer code developed for 


this numerical integration is contained in Appendix C. 


3:5 


Hil, ROADWAY INPUTS 


This chapter explains the road inputs that were used to drive 
the computer model. The roadway inputs are defined in terms of the 
vertical change upward or downward of the values of the road inputs 
U;, U2, and U3 that are shown on Figure 2-10. These are independent 
vertical displacement inputs to tires for each of the wheels in the 
computer model. As defined in Chapter II, there is one input for the 
front HUMMV wheels, (U;), a second input for the rear _HUMMV 
wheels, (U>), and a third input for the howitzer wheels, (U3). The 
model can be driven over any terrain that can be defined in terms of 
the road inputs U; U2 and U3, The road inputs are functions of time 
which allows the user to control the speed at which the model is 
driven by using the time, velocity, and distance relationships. For 
this design study the road inputs used were a washboard course that 
was like one of the courses used in the actual test at Aberdeen 
Proving Grounds, MD (see Chapter VI), and a predicted "worst case” 


pothole road. 


A. WASHBOARD COURSE 

The washboard course that was used to drive the model was 
constructed to be similar to the Six Inch Washboard Course at APG. 
During the instrumented test at APG there were several data _ runs 
taken on the Six Inch Washboard Course (again, see Chapter VI). 


Comparing the model results for this type of washboard course to the 
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actual instrumented data from the APG test for a similar washboard 
provided a means of validating the output results of the computer 
model. 

To model the APG Six Inch Washboard Course a_ sinusoidal input 
was used with an amplitude of three inches and a distance from peak 
to peak of 72 inches. However, because the distance between the 
front HUMMVYV wheels and both the rear HUMMV wheels and also the 
howitzer wheels does not correspond to even multiples of the 72 inch 
peak to peak distance, a lag had to be built into the input sinusoid for 
both the rear HUMMV wheels and the howitzer wheels. This lag was 
necessary, to get the correct timing of the bump wheel inputs for all 
the wheels in the system. 

Additionally, in order to prevent difficulties in the integration 
routine that would occur if the wheel inputs were all turned on at 
the start of the integration run, the rear HUMMV and howitzer wheel 
inputs were not turned on initially. These inputs were turned on at 
the time when the respective wheels would start hitting the 
washboard course bumps if the HUMMV/howitzer system entered 
the course from level terrain. In other words, the front HUMMV 
wheels entered the washboard and _ started taking road input, U), 
while the rear HUMMV wheels and the howitzer wheels remained on 
level terrain which provided zero input to them. Then the rear 
HUMMV wheels entered the washboard and _ started taking road 
input, Uz, while the front HUMMV wheels continued to have the 


input, Uj, and the howitzer wheels continued to have zero input. 
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Finally, the howitzer wheelsentered the washboard and at this point 
all three wheels had their separate inputs U;, U2, and U3. 

This sinusoidal input, for which the wavelength was defined in 
terms of horizontal distance between peaks, was converted into a 
function of time by using time, velocity, and distance relationships. 
To calculate the time between’ successive peaks for a_ specified 
velocity the inputs U ,, U2, and U3 were then defined in terms of time 
which allowed the computer model to be "driven" at any desired 
Speed. Making the input a function of time was required for using 
the time based ODE 45 integration routine of Matlab to solve the 
equations of motion. Figure 3-1 depicts the wheel inputs U;, U2 and 


driven at 5.0 MPH. 





U3 for the washboard course being 





Washboard Wheel Inputs 


Road Ht (in) 


Time (sec) 


Front HUMMYV Wheel Inout, U1 
wines Rear HUMIMV Whee! Inout, U2 
Solid Howitzer Wheel Input, US 


Figure 3-1 Six Inch Washboard Course Wheel Inputs 
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B. "WORST CASE" POTHOLE COURSE 

The model was also "driven" over a trench type pothole to 
Simulate some of the most severe road conditions that the howitzer 
would be subjected to during tactical operations. The model results 
for A-Frame force from this terrain were used to predict if the A- 
Frame and T-Bar were strong enough to handle this “worst case" 
terrain. The pothole course road inputs, U;, Uy, and U3, are 
determined by a separate computer program (see Appendix D)_ that 
computes the changes in axle height as each of the HUMMV wheels 
and the howitzer wheels roll through the pothole. The road inputs 
are based on the path that would be traced out by the center of a 17 
inch radius wheel rolling through the _ specified pothole. The 
computer program that was developed to determine the road inputs 
allows the user to specify the desired dimension of the pothole for 
driving the computer model. For this design study the "worst case" 
pothole was determined to be 20 inches deep and 30 inches across 
with 85 degree slopes on each edge. This size was selected because it 
was judged to be the largest pothole that the HUMMV driver would 
traverse without stopping the HUMMV. Thus, the "worst case” 
pothole dimensions were selected based on what was judged to be 
the most severe terrain the HUMMV/howitzer would — encounter 
while traveling during tactical operations. The "worst case" pothole 
selected was slightly larger than the dimensions of the pothole that 
was used in an informal test conducted at RIA in June 1992. Figure 
3-2 shows the pothole profile, and the axle path traced for the 


howitzer or HUMMV wheel rolling through this “worst case" situation. 
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Figure 3-2 "Worst Case" Pothole Profile and Axle Trace 
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TV. MODEL RESULTS 


This chapter contains the results of driving the computer model 
over the terrain that was explained in Chapter III. The results from 
driving the model over the washboard course show close agreement 
actual results obtained from the APG test. A detailed comparison of 
the model results to the actual results is contained in Chapter VI. 
The model's results from the "worst case" pothole were used _ to 
predict the maximum value of A-Frame force for the howitzer. From 
this predicted maximum force the stress in the A-Frame and T-Bar 
was calculated and compared to the yield stress for the material (see 
Chapter V). 

The predicted A-Frame force is calculated by multiplying’ the 
combined spring stiffness of the A-Frame/T-bar times the change in 
length of A-Frame. The change in length of the A-Frame is a 
function of the difference in the angular displacements of the gun 
tube mass and the howitzer trail/saddle mass at a given time. The 
angular displacements of the gun tube and the trail were determined 
from numerically integrating the equations of motion. 

In addition to determining results in terms of A-Frame force, 
the system's mass and stiffness matrices (see Appendix B) were used 
to determine natural frequencies and mode shapes of the system. 
These natural frequencies and mode shapes were used to_ gain 
physical insight into the system's performance. The natural 
frequencies and mode shapes for the system are shown in Appendix 


ie 
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The Matlab "FFT" function was used to perform the frequency 
analysis on the A-Frame force results from the washboard course to 
determine the frequencies at which the A-Frame force occurred in 
the model results. This frequency analysis a used 1250 Hz data 
Sampling rate. This sampling rate yielded a time increment, dt, that 
equals 1/1250 or the inverse of the sampling rate. In order to have 
a frequency increment that was sufficiently small, 8192 data _ points 
were used for the "FFT". The record length, T, equals dt times the 
number of points or 8192/1250. The frequency increment, df, was 
the inverse of the record length or 1250/8192. 

A. WASHBOARD COURSE RESULTS 

The time histories of the computer model's predictions for A- 
Frame force from the Six Inch Washboard Course are periodic with a 
period that corresponds to the period of the washboard input. But 
unlike the sinusoidal washboard input the A-Frame force output 
predictions are not sinusoidal. These non-sinusoidal A-Frame _ force 
Outputs most likely results from the effects of the non-linear shock 
absorber damping that was incorporated into the model (see Chapter 
Pe 

The frequency analysis of the computer model's predictions for 
A-Frame force shows that as expected the largest component . of 
output force occurs at the same frequency as the washboard input 
frequency. In all cases, there is also a_ significant but lesser force 
component that occurs at a _ frequency that is twice the input 
frequency. 

The computer model was driven over the six inch washboard 


course at speeds of 3.75, 5.0 MPH, and 5.25 MPH. These _ speeds 
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correspond to the speeds for which actual data was taken during the 
APG test. The A-Frame force levels were relatively low at 3.75 MPH. 
The force levels in the A-Frame reached a maximum at 5.25 MPH 
because this speed for the washboard course corresponds to a 
natural frequency of the system. Driving the model at speeds higher 
than 5.25 MPH resulted in decreasing A-Frame force levels. The 
time histories and frequency analysis of the model's predicted 
results for A-Frame force are shown in Figures 4-1 through 4-6. 
I. Results at 3.75 MPH 
The time domain model results at 3.75 MPH are shown in 
Figure 4-1. The predicted A-Frame forces are approximately 1500 
Ibs in tension and 1200 lbs in compression. Figure 4-2 shows the 
frequency analysis at this speed. As expected, there is a large force 
component at the input of 3.75 MPH or a frequency of 0.92 Hz and 
there is smaller force spike at twice this frequency, 1.84 Hz. 
2. Results at 5.0 MPH 
The time domain results are shown in Figure 4-3. The A- 
Frame force levels reach 2800 lbs in tension and 4800 Ibs in 
compression. Figure 4-4 shows the frequency analysis results. Again, 
a large force component at the input, 5.0 Hz, and another force spike 
at twice this frequency, 2.42 Hz. 
3. Results at 5.25 MPH 
The time domain results are shown in Figure 4-5. The A- 
Frame force levels reach 2900 lbs in tension and 4700 lbs in 
compression. Figure 4-6 shows the frequency analysis results. A 
large force component exists at the input frequency, 1.28 Hz, and 


smaller force component exists at twice this frequency, 2.56 Hz. 
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Figure 4-1 Model Results for A-Frame Force at 3.75 MPH (Time 
History) 
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FFT at 3.75 MPH 
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Figure 4-2 Model Results for A-Frame Force at 3.75 MPH (FFT) 
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A-Frame Force 5.0 MPH 
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Figure 4-3 Model Results for A-Frame Force at 5.0 MPH (Time 
History) 
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FFT at 5.0 MPH 
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Figure 4-4 Model Results for A-Frame Force at 5.0 MPH (FFT) 
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A-Frame Force 5.25 MPH 
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Figure 4-5 Model Results for A-Frame Force at 5.25 MPH (Time 
History) 
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FFT at 5.25 MPH 
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Figure 4-6 Model Results for A-Frame Force at 5.25 MPH (FFT) 
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B. NATURAL FREQUENCIES AND MODE SHAPES 

The natural. frequencies and mode shapes of the system were 
determined by using the mass and_ stiffness matrices that were 
developed in Chapter II. The eigenvalues and eigenvectors of the 
system are determined from the matrix that results from multiplying 
the inverse of the mass matrix, [M-!], by the stiffness matrix, [K]. 
The system's natural frequencies are the square root of _ the 
eigenvalues. The system's mode shapes are the eigenvectors that 
correspond to these natural frequencies. Analysis of these natural 
frequencies and mode shapes shows that the howitzer's bounce mode 
occurs at a natural frequency of 1.28 Hertz which corresponds to a 
speed of 5.25 MPH for the howitzer traveling over the washboard. 
By providing an excitation at a natural frequency of the system (1.28 
Hz) the maximum A-Frame force that could possibly be developed 


for this particular washboard input resulted. 


C. "WORST CASE" POTHOLE COURSE RESULTS 

The model results from the “worst case" pothole that was 
denied in Chapter IIT show that maximum A-Frame force occurs 
Ibs at 15 MPH. At speeds slower or faster than this speed the model 
results showed lower levels of A-Frame force for this pothole. Figure 
4-7 shows the A-Frame force for the “worst case" pothole at 15 MPH. 
There are spikes in the A-Frame force when the wheels of the 
system hit the pothole. On Figure 4-7, at 0.1 seconds the front 
HUMMV wheels hit, at 0.6 seconds the rear HUMMV wheels hit, and 


at 1.2 seconds the howitzer wheels hit the _ pothole. When the 
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howitzer wheels hit at 15 MPH the maximum A-Frame force is 8,500 


Ibs. 


A-Frame Force 





Figure 4-7 Model Results of A-Frame Force for "Worst Case" Pothole 
at 15 MPH 
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V. STRUCTURAL ANALYSIS 


This chapter explains the structural analysis that was 
performed on the howitzer in order to determine both the _ force 
levels in the A-Frame and the stress on the howitzer's components. 
To determine the actual force levels in the A-Frame, the A-Frame 
was instrumented with strain gages as shown in Figure 5-1. This 
instrumented A-Frame was’ then’ calibrated to determine a 
calibration factor which would correlate strain gage readings to A- 
Frame force levels. This calibrated A-Frame was used during the 
dynamic test (see Chapter VI) at Aberdeen Proving Grounds (APG) to 


measure the actual force levels in the A-Frame. 


A. A-FRAME CALIBRATION 

The primary purpose of calibrating the A-Frame was to develop 
a reliable means to correlate dynamic strain gage measurements on 
the legs to the amount of force that causes the strain. Additionally, 
the strain measurements obtained during calibration were used to 
determine the value of the equivalent linear spring stiffness of the 
A-Frame for the computer model. During the calibration extra strain 
gages were attached to help identify levels of stress at key locations 
other than the locations needed for the force versus _ strain 
calibration. These extra gages provided information that was used to 


form conclusions about the stress in the A-Frame, discussed later. 
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NOTE: There were ten single strain gages, G1 through G10, 
located as indicated on the diagram. For the rosette, G11 
is the transverse leg, G12 is the 45 degree leg and G13 is 
the axial leg. 





Figure 5-1 Location of A-Frame Strain Gages 
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The A-Frame was loaded in both a _ vertical and lateral load 
configuration to determine if there were any significantly different 
effects between the two loading configurations. Figures 5-2 & 5-3 
show the test setup for both these configurations. 

i A-Frame’ Instrumentation 

To conduct the calibration, the A-Frame was _ instrumented 
with 13 strain gages (Figure 5-1). The data taken from the pair of 
single strain gages located on the bottom of each leg was used for the 
development of the calibration curve. The pair of single strain gages 
on the top of each leg were used to verify that loads in the legs were 
mainly axial loads as opposed to bending loads. The single gage on 
the A-frame cross support brace and the single gage on the attaching 
Clamp for this brace were used to identify stress levels in_ these 
members. The latter two strain gages were used because there were 
some known cases where the cross support brace clamps _ had 
cracked, and it was unknown whether these cracks were the result of 
high load levels or caused by manufacturing defects. 

pap Support Equipment 

The calibration was conducted using the lab facilities of the 
Department of Mechanical Engineering at the Naval Postgraduate 
School. To hold the A-Frame during the calibration a_ test support 
structure was built. The main objective of this support stand was to 
closely match the mechanical effects that occur in the A-Frame when 
it is loaded while connected to the T-Bar. The test support used was 
constructed out of heavy steel and used actual eyebolts from a M119 


Howitzer to connect the bottom of the A-Frame legs to the test 
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Figure 5-2 Vertical Configuration for A-Frame Test 
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Figure 5-3 Lateral Configuration for A-Frame Test 
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support (Figure 5-4). A clamping device to attach the apex of the A- 
Frame to the load cell was also designed which simulated the loading 
effects of the actual clamping device on the howitzer (Figure 5-5). 
The test support structure achieved the goal of simulating the A- 
Frame being attached to the T-Bar of the howitzer. 
3. Load Cells 

The vertical loads were applied using a Riehe Testing 
machine as shown in Figure 5-2. This loading device was used to 
apply vertical loads up to 12,000 Ibs. The lateral loads were applied 
using a _turn-buckle that was connected to a 5,000 pound rated 
capacity dynamometer. This dynamometer and_ turn-buckle system 
is shown in Figure 5-3. 

4, Calibration Procedures 

Two data runs were recorded for both the vertical and 
lateral configurations. A data run consisted of loading the A-Frame 
in increments to the maximum load for that particular data run and 
then unloading it in the same manner. All strain gages were zeroed 
before the start of each run. At each load increment, strain gage 
readings were recorded for all of the strain’ gages. Tables F-1 
through F-4 contain the strain gage readings and load data obtained 


from these calibration procedures. 
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Figure 5-4 A-Frame Connection to Test Mount 
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Figure 5-5 A-Frame Apex Clamping Device 
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5. Calibration Factor 
From the vertical test strain gage data in Appendix F _ the 
A-Frame vertical force versus’ strain calibration factor was 
determined. Since the A-Frame vertical load test was a static. test 
the vertical force loading at the apex of the A-Frame must be 
balanced by the sum of the vertical reaction forces in the legs of the 


A-Frame (Figure 5-6). 


Strain | Geeaini 
Gage Pair #1 Gage Pair # 4 


a 
/ 
28 


\ 
Ry Nv2 28° 
where: 
F = the applied vertical load at the apex of the A-frame 


Ri & Rio = the horizontal reaction forces where the A-frame 
connects to the T-BAR 


Ry; &R» = the vertical reaction forces where the A-frame 
connects to the T-BAR 





Figure 5-6 Balance of Forces on the A-Frame 
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The reaction forces (Ry; &Ry72) at the bottom of each leg of 
the A-Frame were determined by using the strain gage readings 
from Gage Pair 1 (gages 1 & 2) and Gage Pair 4 (gages 9 & 10) that 
are identified in Figure 5-6. The average value of the strain for 
either of the strain gage pairs was multiplied by the Modulus of 
Elasticity of the A-Frame (10 x 10® psi) which yields the axial stress 
in the A-Frame leg. This axial stress value was then multiplied by 
the cross. sectional area of the A-Frame legs (0.6107 in-) to 
determine the axial force. To calculate the vertical component of this 
axial force the above value for axial force was multiplied times the 
cosine of 28°. A sample calculation is shown below using the 10,000 
Ib load increment of Vertical Test #1 (see Table F-1 for strain gage 


data): 


Strain Gage Pair 1 = Avg Strain = 988 10° 
.. Stress = Strain x Modulus of Elasticity 
= (988 x 10°°)(10 = 10°psi) = 9,880 psi 
= Force = Stress x Area = (9,880 psi)(0.6107 in’) = 6,033 Ibs 
. The vertical reaction force: R,, = (6,033 lbs)(cos(28)) = 5,327 Ibs 


Strain Gage Pair 4 = Avg Strain = 1,095 «10° 
. Stress = Strain x Modulus of Elasticity 
= (1,095 x 10° )(10 x 10°psi) = 1,095 psi 
= Force = Stress x Area = (10,950 psi)(0.6107 in’) = 6,687 Ibs 
.. The vertical reaction force: R,, = (6,687 IbsXco¥28)) = 5,904 Ibs 





.y’ Vertical Forces = R,, +R,, = 5,327+ 5,904 = 11,231 Ibs 
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To determine the calibration factor for vertical loading, 
the percentage difference between the known apex load and the sum 
of the calculated vertical force components at the bottom of the A- 
Frame was determined. The average percent difference for all the 
load increments is equal to the calibration factor. This calibration 
factor can then be used to compare measured loads in the A-Frame 
legs to the applied load at the apex of the A-Frame. An_ average 
calibration factor of 0.88 was determined over all of the load 
increments from both vertical tests. A sample calculation of the 
calibration factor for the 10,000 lb force increment of Vertical Test 
#1 (see Table F-1) is shown below: 

Actual Vertical Force = 10,000 Ibs 


Measured Vertical Force = 11,231 lbs 


Pabrinone Eactore Seo 


11,231 


6. A-Frame Stiffness 

The measured strain data was also used to determine an 
effective linear spring stiffness for the A-Frame. To determine this 
stiffness the applied vertical load for a specified load increment was 
divided by the corresponding amount of vertical deflection in the A- 
Frame. The average axial deflection in the A-Frame legs was 
calculated from the average strain at the bottom of each leg. The 
average axial deflection was then multiplied by the cosine of 28°. to 
calculate the amount of vertical deflection. A sample calculation for 
the 10,000 lb load increment of Vertical Test #1 (see Table F-1) is 


shown below: 
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Strain Gage Pair 1 = 998; 


Avg axial strain oe 1042 x10°°; Multiply by leg length 


Strain Gage Pair 4 = 1095 






= Avg axial deflection = (28 in)(1042 x 10°) = 0.029 in 
Convert to vertical deflection = (0.029 in )(cos(28)) = 0.0257 in 


10,000 Ibs 
0.0257 in 





.. A-Frame Stiffness = = 388,400 lbs/ in 





ade Observations 
From analyzing the results of both the vertical and lateral 
A-Frame loadings the following conclusions were reached regarding 


A-Frame loading: 


* The strain measurements in both the cross member brace and 
the attaching clamp for this brace on the A-Frame (strain gages 3 
& 4, in Figure 5-1) were very low regardless of the force at which 
the A-Frame was loaded. This indicates that these members are 
not critical components of the A-Frame in regards to failure. 


* The A-Frame loads are carried mostly through axial loads in the 
A-Frame legs as opposed to bending loads. This means that the 
bottom connection joints where the A-Frame attaches to the T-Bar 
act as pinned joints However, because there was some bending in 
the A-frame legs at lower force levels the bottom joints do not act 
as perfect pin joints. The higher force levels break the static 
friction effects in the joints causing the joints to act more as 
perfect pin joints. 


B. STRESS ANALYSIS 
The stress analysis focused on determining the highest stress 
levels in the T-Bar and the A-Frame. To determine if mechanical 


failure will occur due to exceeding the yield stress, the stress at the 
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probable points of highest stress in the A-Frame and T-Bar were 
compared to the yield stress values for the respective materials. 
1. T-Bar Stress 

The T-Bar and A-Frame assembly are shown below in 
Figure 5-7. Point A on this figure is identified as the point of highest 
stress. This is where the geometry of the T-Bar assembly changes 
and the stronger built up portion of the T-Bar ends. To simplify the 
Stress analysis the T-Bar is modeled as a cantilevered beam. Figure 
5-8 below depicts the cantilevered beam modeling of the T-Bar as 
viewed from the end of the T-Bar and the rear of the howitzer. The 
point of highest stress was determined to be at the corner of the T- 
Bar at the base of the modeled cantilevered beam. One of these 


corners is identified as Point A in Figure 5-8. 


{ao eae A-Frame 


Transverse | 


Transverse 


Hole 
Hole 
High Stress Lo High Stress 
Location Location 


(Point A) 


=—s ZN 


8 inches ‘8 inches’ 





Figure 5-7 T-Bar and A-Frame Assembly 
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View from Rear of Howitzer View from end of T-Bar 


Point A 





Figure 5-8 T-Bar Cantilevered Beam Model 


To calculate the stress at Point <A, the normal _ stress 
components that result when the axial force in the A-Frame leg, P, 
was resolved into components in the x, y, and z directions as shown 
on Figure 5-8. The stress at Point A is the sum of these -normal stress 
components using the principle of superposition of stress. To find 
the normal stress components first the A-Frame leg force, P, was 
projected into a force component along the axis of the T-Bar in the y 
direction, Py = sin(28)P, and a second force component in the xz 
plane, Pxz= cos(28)P. The P,zcomponent was further broken up into 
components in the x and z directions: Px, = cos(27)Pxz and Poe— 
sin(27)P,7. To calculate the normal stress that was caused by _ these 
force components Px, Py, and Pz well known strength of materials 
formulas for stress in a beam are applied. Therefore, applying 
these formulas the expression for the maximum stress in the T-Bar 


1S: 
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ECE) oo 


A I Nt 


X Zz 


where: 


Px, Py, and P,: are defined in the above paragraph. 


L=8 inches: The distance from the base of the cantilever to the 
point of application. 


A=.9240 in*: The cross sectional area of the T-Bar. 
Ix=.8962 in*: The bending moment of inertia about the x axis. 


I1,=.5979 in*: The bending moment of inertia about the z axis. 


x =1.300 in: The distance in the x direction from the geometric 
center of the T-Bar to the outer edge of the T-Bar. 


z =1.005 in: The distance in the z direction from the geometric 
center of the T-Bar to the outer edge of the T-Bar. 


Pah A-Frame Stress 
The probable point of highest stress in the A-Frame. will 


occur around the 0.375 inch diameter transverse circular hole that 


positions the cross member attaching clamps. This hole is identified 
in Figure 5-7 above. At the location of this hole a_ stress 
concentration will exist. Figure 5-9 below contains the empirical 


formula from page 601 of [Ref. 4] that was used to determine the 
stress concentration in a pipe section with a transverse hole. 

From the above formula the stress concentration factor for 
axial loads owas determined to be = 3.73. Additional stress 
concentrations due to the effects of bending and torsional loads were 


neglected. The justification for neglecting the effects of bending 
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Radial hole in a hollow or sold | 18a Ejastic stress, axial tension | = ia a where £ = Ky + K,(=) + K(=) + Kf 


| | and where for d/D ¢ 0.9 and 2r/D ¢ 0.45 
| 
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3.000 
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1 Ky = —6.793 + 1.1334/D — 0.126(d/D)? 
38.382 — 7.242d/D + 6.495(d/D)? 


Ky = —44.576 — 7.428d/D + $8.656(d/D)? 





Figure 5-9 Stress Concentration Due to a Circular Hole in a Pipe 
Section 


loads was that both the A-Frame calibration results (see Appendix F) 
and the results of the instrumented test (see Chapter VI) show that 
the joint where the A-Frame connects to the T-Bar acts as a nearly 
perfect pin joint at higher load levels. Any stress concentration 
reduction due to the spring pin that is driven into the hole were 
neglected because this stress reduction was assumed to be_ small. 
Therefore, maximum total stress in the A-Frame was based on only 
the axial stress component in the legs and was determined using the 


following formula: 


PK 
oe. (Eqn 5-2) 
where: 
P = the axial load in the A-Frame leg. 


K = 3.73: the stress concentration factor for axial loading that is 
caused by the hole in the A-Frame leg. 


A = 0.6107 in?: the cross sectional area of the A-Frame leg. 
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Without this hole in the A-Frame the maximum _ allowable 
Stress would be increased by the amount of the stress concentration 
factor, 3.73. This hole appears to only aid in assembling the A-Frame 
and serves no other function. Elimination of this hole would make 
the factor of safety more than three times greater. 

5; Worst Case Stress and Factors of Safety 

Using Eqn 5-1 and Eqn 5-2 above, the maximum stress was 
calculated for the T-Bar and the A-Frame. The maximum predicted 
A-Frame force of 8,500 ibs from driving the computer over the 
“worst case" pothole described in Chapter II was divided by two to 
calculate the force in each leg of the A-Frame. The resulting axial 
force, P, in each leg of the A-Frame equals 4,250 lbs. From Eqn 5-1 
the resulting stress in the T-Bar is 63,900 psi. From Eqn 5-2 the A- 
Frame stress equals 26,000 psi. 

The above stress levels are used to calculate the predicted 
factors of safety for the T-Bar and the A-Frame. The factor of safety 
is calculated by dividing the yield stress of the component by the 
predicted maximum load stress. The yield stress for the T-Bar and 
A-Frame are 116,000 psi and 50,000 psi, respectively. Based on the 
maximum predicted stress levels given above, the factor of safety is 
1.81 for the T-Bar and 1.92 for the A-Frame. These factors of safety 
show that both the T-Bar and the A-Frame are strong enough to 


withstand the stress from the worst case predicted loads. 
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VI. INSTRUMENTED TEST 


An instrumented dynamic test was conducted on the M119 
howitzer at Aberdeen Proving Grounds (APG) in April of 1993. The 
howitzer was instrumented as described below in paragraph A. The 
dynamic testing was conducted by driving the howitzer over several 
of the APG courses at various speeds. The calibrated A-Frame 
described in Chapter V was used to measure the force levels in the 
A-Frame during the test. In most respects, the test results for A- 
Frame force are similar to the results obtained from the computer 
model. However, because certain physical effects were not included 
in the computer model, there are some differences between the 
results of the actual test and the results from the computer model. 
The results of this instrumented test confirmed that the M119 
Howitzer could be safely transported in the firing position with the 


A-Frame attached. 


A. INSTRUMENTATION 
For this test the howitzer was instrumented with 20 © strain 
gages, two tri-axial accelerometers, and two _ extensiometers. This 
instrumentation is explained in the following paragraphs. 
1. Accelerometers 
The howitzer was instrumented as shown in Figure 6-1. A 
tri-axial accelerometer was placed on the T-Bar to measure the 


vertical, longitudinal, and transverse accelerations of the howitzer 
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Figure 6-1 Howitzer Instrumentation for APG Test 








trail mass. This accelerometer was not placed at the center of 
gravity of the trail mass because of difficulties in attaching an 
accelerometer at this _ point. The location chosen on the T-Bar 
provided accelerations on the same order of magnitude and 
representative of the accelerations at the center of gravity of the 
trail. The other tri-axial accelerometer was attached on the elevating 
mass near its center of gravity. This location is shown in Figure 6-2. 
Je Single Strain Gages 

The two pair of single strain gages on the bottom of each leg 
(gage pair 1 & 4 in Figure 5-6) of the previously calibrated A-Frame 
were used to measure the strain in the A-Frame legs. The calibration 
factor determined in Chapter V was used to correlate measured 
Strains in the A-Frame legs to force levels in the A-Frame. Single 
Strain gages were also attached on the top and the bottom of the T- 
Bar on both sides of the howitzer. The location of these gages is 
shown in Figure 6-3. Because the highest levels of stress on the T- 
Bar are at the corners, it was not possible to place the gages at these 
points. The gages were placed as close to the high stress points as 
possible. Even though these gages could not be used to determine 
the maximum stress on the T-Bar, they could be used to roughly 
judge the level of stress. 

3. Other Gages 
Strain rosettes were placed on the bottom of the A-Frame 


legs and on the bottom howitzer trail on each side of the howitzer. 
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Figure 6-2 











Location of T-Bar Strain Gages 


Figure 6-3 
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An extensiometer was used on each shock absorber to determine if 
the shocks were extending far enough so that the stiff rubber bump 


stops would be activated. 


B. TEST COURSES 
Several of the standard APG test courses were used for the 
instrumented test. These courses are described below. 
i Belgian Block Course 
This course is designed using "Belgian Blocks" to produce a 
road that has high frequency random bumps. The magnitude of each 
bump is small. This course is shown in Figure 6-4. During the test 


this course was traversed at speeds ranging from 10 to 25 MPH. 


2° BEO OF SANO 8& CEMENT 6" CONCRETE 


— GRANITE 3LOCKS IN CONCRETE 6 CONCRETE 


LENGTH OF “COURSE S.94C0F] 





Figure 6-4 APG Belgian Block Course 


Zz. Six Inch Washboard Course 
This course has six inch bumps that are evenly spaced six 


feet or 72 inches apart as shown in Figure 6-5. This course produces 
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Significant heave and pitch effects on wheeled vehicles even 
relatively slow speeds. During the test this course was traversed 


speeds ranging from three to six MPH. 


LONGITUDINAL SECTION 


LENGTH OF COURSE, 798 FT. 


Figure 6-5 APG Six Inch Washboard Course 


3. Cross Country Course Number 1 


at 


at 





This course is similar to a relatively smooth but winding 


trail road. This course was traversed at approximately 25 MPH. 


4. Cross Country Course Number 3 


This 1s similar to a rough trail road which has random large, 


medium, and small sized bumps. This course was _ traversed 


speeds ranging from 5 to 15 MPH. 
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C. TEST PROCEDURES 

Each of the above courses were traversed by the howitzer being 
towed by the HUMMV at speeds within the ranges listed. Data was 
collected on all the channels described above, using a sampling rate 
of 2,500 Hertz. Because of memory limitations in the data collection 
equipment only approximately 30 seconds worth of continuous’ data 
could be stored. For ‘each course several data runs were completed 
using various speeds. For instance, on the Six Inch Washboard 
Course, six data runs were taken at speeds ranging from three to six 
MPH. 

The speeds for each data run on the Six Inch Washboard Course 
and Cross Country Course Number 3 were measured using the 
HUMMV's speedometer, because these courses were too rough to use 
more accurate devices. For the data runs, on these courses, the 
driver tried to hold the HUMMV at the specified speed while himself 
being bounced around. Because of the extremely rough road 
conditions and inaccuracies in the speedometer, the actual speeds 
ended up being different than the speeds specified in the test plan. 
Consequently, the only accurate method of determining the actual 
speed during a data run on these courses is from the video tape 


taken during the run. 


D. TEST RESULTS AND DISCUSSION 
The Belgian Block and Cross Country Course Number 1 showed 
no significant levels of strain or acceleration even at speeds up 25 


MPH. These tests simply proved that the howitzer could traverse 
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these courses without causing damage to the howitzer. The Six Inch 
Washboard Course and Cross Country Course Number 3, however, did 
produce significant levels of strain and acceleration on the howitzer 
components. 

The test results for the Six Inch Washboard Course’ were 
analyzed in both the time and frequency domains. Data runs at 
speeds of 3.75, 5.0 and 5.25 MPH were analyzed for the Six Inch 
Washboard Course and compared to the predicted results from the 
computer model. The most severe data run from Cross Country 
Course Number 3 at 15 MPH was analyzed in the time domain. The 
procedures used and the analysis of these results are as follows: 

ie Time History Results 

The time history for vertical "in plane" A-Frame force in 
each case was produced by first summing the “in sane axial forces 
from both of the A-Frame legs. The same calculation procedures as 
in Chapter V were used for determination of the axial forces in the 
legs from the strain gage data for strain gage pairs 1 and 4 (see 
Figure 5-6). The vertical force was determined by multiplying the 
axial force times the cosine of 28°. This total vertical force in the A- 
Frame legs was multiplied by the calibration factor of 0.88 that was 
determined in Chapter V. The A-Frame lateral forces were 
calculated by determining the difference in the horizontal force 
components of the axial force in each A-Frame leg. The accelerations 
were plotted directly from the measured data. 

The force and acceleration time histories for the washboard 


course (Figures 6-6 through 6-11) show that high frequency spikes 
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occur on the washboard course. These spikes are caused the braking 
System being activated as the howitzer went over the bumps for all 
speeds. Also, for the 5.25 MPH data run the spade hit the tops of the 
bumps. These effects are not included in the computer model and 
result in differences between the model results and the actual 
results, particularly, for the 5.25 MPH data because the spade was 
hitting. 

Time histories for A-Frame force for speeds 3.75, 5.0 & 5.25 
MPH from the actual test (Figures 6-6, 6-8, & 6-10) are similar to the 
time histories predicted by the computer model for the same _ speeds 
(Figures 4-1, 4-3 & 4-5). The high frequency force spikes do _ not 
appear in the model's predicted results, however, the levels of the 
basic force response in the actual results closely compare to the force 
response of the model. For example, the actual results or 5.0 MPH 
(Figure 6-8) has spikes that reach a maximum force of 5,500 Ibs. But 
the basic force response on average reaches approximately 2,400 lbs 
in tension and 4,300 Ibs in compression. The corresponding model 
prediction (Figure 4-3) shows forces of 2,800 lbs in tension and 
4,800 lbs in compression. A _ similar comparison can be made at the 
3.75 and 5.25 MPH speeds. Although ne actual results are not 
exactly the same as the model predictions, they are close enough to 
conclude that the model is reasonably representative of the actual 
results. 

The A-Frame vertical force results for the most severe 
portion of the Cross Country Course Number 3 are shown in Figure 6- 


12. These results show that on a rough trail road the A-Frame force 
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is not as great as the force that was caused by the Six Inch 


Washboard Course at a speed of 5.25 MPH. 


Vertical Force from APG Test (3.75 MPH) 
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Figure 6-6 A-Frame Force at 3.75 MPH on the Washboard Course 
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Vertical Trail Accel from APG Test (3.75 MPH) 
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Figure 6-7 Trail Vertical Acceleration at 3.75 MPH on_ the 
Washboard Course 
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A-Frame Vertical Force from APG Test (5.0 MPH) 
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Figure 6-8 A-Frame Force at 5.0 MPH on the Washboard Course 
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Vertical Trail Accel from APG Test (5.0 MPH) 
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Figure 6-9 Trail Vertical Acceleration at 5.0 MPH on the Washboard 
Course 
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A-Frame Vertical Force from APG Test (5.25 MPH) 
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Figure 6-10 A-Frame Force at 5.25 MPH on the Washboard Course 
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Figure 6-11 Trail 
Washboard Course 


Vertical Trail Accel] from APG Test (5.25 MPH) 
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Vertical Acceleration at 5.25 MPH 
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Vertical Force for Cross Coutry #3 (15 MPH) 
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Figure 6-12 A-Frame Force at 15 MPH on Cross Country Course 
Number 3 
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rae Frequency Domain 

The frequency domain analysis for the actual A-Frame 
force and trail acceleration results are shown in Figures 6-13 through 
6-18. The same procedures using the "FFT" function in Matlab as 
were used in Chapter IV for the model results were used to perform 
the frequency analysis on the actual _ results. The time increment 
used was 1/1250 seconds and was determined by using every other 
data point of the original data that was taken at 2500 Hz. 8192 
points were used to achieve an adequate frequency discretization. 
The record length was 8192/1250 and _ frequency increment’ was 
1250/8192. The FFT analysis of both A-Frame force and_ trail 
vertical acceleration for each data run clearly shows two frequency 
spikes. The lowest frequency and largest spike corresponds to the 
forcing frequency or road input to the howitzer tires. The second 
largest spike corresponds to twice the forcing frequency. 

The frequency analysis for actual A-Frame force for speeds 
of 3.75, 5.0, & 5.25 MPH (Figures 6-13, 6-15, & 6-17) compares 
reasonably well to the corresponding frequency analysis from_ the 
model results (Figures 4-2, 4-4, & 4-6). For example, the FFT results 
from the actual data at 5.25 MPH (Figure 6-17) shows a large force 
component 1180 lbs at the road input frequency (1.28 Hz) and 
second large force component of 575 lbs at double the road input 
(2.56 Hz). Where as the FFT of the model for the corresponding 
predicted A-frame force (Figure 4-6) shows a large force component 


of 1375 lbs at the road input frequency (1.28 Hz) and a second large 
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force component of 500 Ibs at twice the input frequency (2.56 Hz). 
For this case, the total amount of force from the 1.28 and 2.56 Hz 
spikes for the actual results and the model results are approximately 
the same (1765 Ibs versus 1875 lbs). Although FFT results for the 
actual data and the model predictions are not exactly the same, they 


are close enough further to validate the model's predicted results. 


Vertical Force FFT from APG Test (3.75 MPH) 
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Figure 6-13. FFT of A-Frame Force at 3.75 MPH on the Washboard 
Course 
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Vertical Trail Accel FFT from APG Test (3.75 MPH) 
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Figure 6-14 FFT of Trail Verticai Acceleration at 3.75 MPH on the 
Washboard Course 
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Vertical Force FFT from APG Test (5.0 MPH) 
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Figure 6-15 FFI of A-Frame Force at 5.00 MPH on the Washboard 
Course 
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Vertical Trail Accel FFT from APG Test (5.0 MPH) 
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Figure 6-16 FFT of Trail Vertical Acceleration at 5.0 MPH on _ the 
Washboard Course 
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Vertical Force FFT from APG Test(5.25 MPH) 
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Figure 6-17 FFT of A-Frame Force at 5.25 MPH on the Washboard 
Course 
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Vertical Trail Accel FFT from APG Test (5.25 MPH) 
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Figure 6-18 FFT of Trail Vertical Acceleration at 5.25 MPH on _ the 
Washboard Course 
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3. Other Results 

The strain reading from the strain’ rosettes that were 
attached to the bottom the howitzer's. trail showed low _ strain 
readings for all courses. These low strain readings clearly showed 
that transporting the howitzer in the firing position with the A- 
Frame attached produces insignificant levels of stress in the howitzer 
trail. 

Lateral forces in the A-Frame were also determined from 
the actual test data. Figures 6-19 through 6-21 show the lateral 
force results for 3.75, 5.0 & 5.25 MPH, respectively. The maximum 
lateral forces were less than ten percent of the maximum vertical 
forces for all runs on the washboard course. For the rough terrain 
Cross Country Course Number Three the maximum lateral force 
(Figure 6-22) was twenty percent of the maximum vertical force. 
These relatively low values of lateral force help to show the validity 
of the modeling technique used, since the computer model _ only 
predicts the vertical forces in the A-Frame and assumes that the 


lateral forces would be comparatively low. 
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A-Frame Lateral Force from APG Test (3.75 MPH) 
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Figure 6-19 Lateral Force from the APG Washboard Course at 3.75 
M PH 
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A-Frame Lateral Force from APG Test(5.0 MPH) 
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Figure 6-20 Lateral Force from the APG Washboard Course at 5.0 
MPH 
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A-Frame Lateral Force from APG Test (5.25 MPH) 
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Figure 6-21 Lateral Force from the APG Washboard Course at 5.25 
MPH 
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Figure 6-22 Lateral Force from Cross Country Course Number Three 
at 15 MPH 
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VII. CONCLUSIONS/IRECOMMENDATIONS 
A. CONCLUSIONS 

The computer model developed, produced output results that are 
reasonably representative of the actual results for the washboard 
type road input used at Aberdeen Proving Grounds. Therefore, the 
model is validated and can be used to predict "worst case" loading on 
the howitzer's critical parts. These "worst case" predictions for a 
trench type pothole show the M119 Howitzer can be safely towed in 
the firing position with the A-Frame _ attached. 

As an additional result of this study, we believe that the T-Bar is 
the most critical component on the howitzer in regards to mechanical 
failure when towing in the firing position. The predicted factor of 
safety for the T-bar based on the “worst case" pothole terrain is 1.81. 
Due to the high cost of repair, making the T-Bar out of thicker steel 
would provide an increased factor of safely. 

The predicted factor of safety for the A-Frame, from the same 
"worst case" terrain, is 1.92. This indicates that the A-Frame is 
strong enough to withstand the worst case expected loads. However, 
the A-Frame is made out of aluminum and aluminum does not have 
a fatigue limit. Thus, eventually the A-Frame is expected to fatigue. 
The A-Frame's fatigue life could be greatly extended by _ eliminating 


the holes in its legs. 
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The strain reading taken from the trail rosettes during the 
instrumented test at APG show tthat the stress in the trail is 


insignificant during towing in the A-Frame position. 


B. RECOMMENDATIONS 

The US Army should remove the towing restrictions and allow 
the howitzer to be towed in the firing configuration with the A- 
Frame attached during tactical operations. 

The T-Bar should be made out of thicker steel to provide added 
assurance that it will not fail. 

The A-Frame should be _ re-configured for manufacture without 


the hole in its legs. 
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APPENDIX A 


SYSTEM PARAMETER VALUES 


The following table defines the HUMMV/Howitzer = system 
parameters using an inches (in), pounds (ibs), and seconds (sec) 
system of units. Gravity (G) is approximated to be equal to 386.4 
in/sec*. The nomenclature for the parameters that are identified in 
this Appendix is the same nomenclature that is used throughout this 
paper. Refer to Figures 2-1 and 2-2 to identify the parameters 
defined below. The values listed in the table incorporate the effects 
of both sides for those components that are identical on the right and 
left side of the HUMMV or Howitzer. For example, the spring 
stiffness for each of the front HUMMV tires, K,, is 1323.0 lbs/in. This 
value is multiplied by two to account for the combined effect of both 
front HUMMV front tires. The value listed in the following table for 


K, is 2646 Ibs/in. 
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TABLE A-1 COMPUTER MODEL PARAMETER VALUES 


How _Determined/Notes 


2646.0 Ibs/in Reference 1/HUMMV front tires stiffness 
1908.0 Ibs/in Reference 1/HUMMV front suspension stiffness 
3544.0 Ibs/in Reference 1/ HUMMV rear tires stiffness 
5040.0 Ibs/in Reference 1/ HUMMV rear suspension stiffness 
ee 2000.0 Ibs/in Estimated /Howitzer tires stiffness 

i erence . Modeled Howitzer vertical 
| K, | 1056.0 Ibs/in nelerene 
ce 200,000.0 Ibs/in | Calculated/Combine A-Frame & T-Bar stiffness | A-Frame | Calculated/Combine A-Frame & T-Bar stiffness | T-Bar stiffness 


Bic=171.6 Ibs-sec/in | Reference 1/Bic is for compression & BiT 
Bit=308.4 Ibs-sec/in | for tension. Values are 20% more than Ref 1 
| Be ee 2 Ibs-sec/in | Reference 1/ B2c is for compression & B2tT 
; i for tension. Values are 20% more than Ref 1 
a Calculated from RIA data/ B3c is for 
compression & B3T for tension. 
2 250,000 Ibs-in-sec/rad | Estimated/ Torsional Damping in Sai clutch 
Re 


1.294 | 1.294 om 


a be si inclicies meee at HUMMV rear 
| 1.294 bm wheels,shocks and suspension springs 
= 21.22 lbm Reference 1/ HUMMV Body 


1.035 Ibm re e 2/Includes ‘el wheels, 


he r a 
a 9.834 ibm none 2/ Howitzer trai & = mass 
ts. 5.952 Ibm Reference 2/ Howitzer elevating mass 


52,680.0 Ibs-in-sec? | Reference 1/ Moment of inertia about HUMMV CG 
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TABLE A-1| COMPUTER MODEL PARAMETER VALUES (CONT) 


oe) 


ymbo Value/units How Determined/Notes 


. Estimated/ Moment of inertia of Trail & saddle 
i 10,000 Ibs-in- secé pean ae 
NTO OGmiboinscecs Estimated/ Moment of inertia of Elevating mass 
: apout (! 
.8 in eference MMV CG to front springs & shocks 
L, 83.8 Ref 1/HU CG to fr & shock 
L, Reference 1/HUMMV CG to rear springs & shocks 
Reference 1/HUMMV CG to hitch connecting point 


130.5 in Measured/ HUMMV front axle to HUMMV rear axle 
294.5 in Measured/ HUMMV front axle to Howitzer axle 


iE BAa0 in Measured/Trunion to where A-Frame 
connects to the T-Bar 
110.3 in Measured/Trunion to hitch connecting point 
ie. 45.0 in Measured/ Trunion to w ere A-Frame 
tt tat ies @ ie oY 3 ° é @,< 
provide by RIA/ Truton 05 of vat mas 
Measured/ Static angle from horizontal t 
Ca 12.0 degrees oa Tubelvnen Aa ae poe " i | 
a 51.0 degrees Measured/ Angle from Gun Tube to A-Frame 
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APPENDIX 8B 


MATRICES FOR EQUATIONS OF MOTION 


This Appendix 


contains 


the 


mass, 


damping, and _ stiffness 


matrices for the equations of motion that were developed in Chapter 


II. The mass system matrix, [M], is as follows: 


Im, 0 0 0 0 0 0 
lo m, 0 0 0 0 0 
lo 


| 0 m; O m3, M3, M;, 
is 0 0 m,, 0 0 0 
10 0 m;;, 9 ms, M5, Ms; 
0 0 m,; 0 M5 Mg 
OO ine 0 eee a 
where: 

m,,=M, m,,=M, 


m,;=M,;L, + M,L, 

m,, =M,L,,, cox @,) 
m,;=M;L, + M,L, 

ms, =M.L,L,, + M,L,L,, 
M;L,, + ML, 

Mm, =I, + ML; + MoLig 
m,, = M,L,,, co%9,) 

M7, = MeL cgli¢ cos(9g) 


M3 = 


Note due _ to 


equal: 


M5; — M,5 M63 
Mes = M56 Mm, 


7 


symmetry of [M] 


m,;—=M;+M,+M, 
m3, =M.L,, + M.L,. 


M44 = 
ms; =I, + MSL, 
m,;, = M,L,,L; cos(0,) 


M, 
t Moles 


Mos =M.L,L,, + M,.L;L,. 


the following matrix 


= N36 


5 M<7 
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m,, = M,L,.,L,,cos(6,) 
m,, =M,L.,L, cos(9,) 
m,, =I, + M Log 


elements are 


M5; ~ M37 


M7, — M67 


The damping matrix for the system is dependent on _ whether 
the shock absorber forces are tensile or compressive at a given time. 
The tensile or compressive state for each of the shock absorbers is 
determined by the relative velocities of the end connection points for 


the shock absorber. The damping matrix, [C(x)], is as follows: 


ion «OO Ce 0 Ce 0 0 
10 2 23 0 C25 0 0 
0 


C3, C35 C33 C3, C35 C36 


where: 

¢,, = B, c,,= —B, ¢,;= BL, 
Co, = B, Co; = ~B, c,, = ~B,L, 
C3, = —B, C3, = —B, c,, = B, + B, + B, 
C,, = —B, c,,=—B,L, +B,L, + BL, 

C35 = B,L,, C,, — —B, C,, = B, 

C,5 = —B,L, C45 = ~B,L,, ¢,,= BL, 
C5, = —BL, c,,= —BL, +B,L, + BL, 

Cs, = —B,L, C..=BL, + BL + BL 

Cs, = BLL, Co, = BL, Ce, = —B3L,, 
C.5= BIL,L, Sayre eee ae JE Co7 = ~B, 
C7, = —B, C7, = B, 


Note due to symmetry of [C(x)] the following elements are equal: 


C3, 7 3 Go = Co; C5, ~ Gs C43 — C34 
Cs = Cos C53 — C35 C54 ~ Cys Ce3 — (36 
Cog — Cae Ces — O56 Cr6 = Coz 
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The system stiffness matrix, [K], is as follows: 


fe 0 Os kn | 

10 =k ke 0 ky. mONOn | 

: kx, ky ks, k,,  k;s k,, 0 | 

ie ke ke k,, O | 

| ks, ks, Ks, ks, = Ks ks, 9 | 

| 0 0 k., ke, Kes Kes Key | 

lo 0 0. 07-8 Seater 

where 

k,,=K, + K, k,,=-K, k,; = K,L, 
k,,=K, + K, k,,= —K, k,, = -K,L, 
k,, = —K, k,, = -K, k,,-K,+K,+K, 
Ke al kK, = oKGL tek eee 

k,, = K,L,, k,,=—-K, k,,=K, + K, 
k,, = -K,L, k,.= -K,L,, ks, = K,L, 

ke sake ee k= —K Dei Eker. 

ks, = —K,L, k,,=K,L,° + K,L,° + KL,’ 

k,, = K{L,L,, k,, = K,L,, k,, = KL, 
k,, = K,L,L, keg =K,L,z + K,L,L,,cos(0,) 

k= = Ke Lacoste) ki. = = Ke ee cosioa) 


ko = Kee costae) 


Note due to symmetry of [K] the following matrix elements are equal: 


Kis Se k,, =k,, k,, =k,, Keen 
k,, = k,. k,, = k,, k,, = k,; k,3 = ky, 
k, = Ky, Kos = Kee ki, = kg, 
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APPENDIX = C 


NUMERICAL INTEGRATION COMPUTER CODE 


This Appendix contains the computer code used to integrate 
the equations of motion using the ODE 45 integration routine in 


Matlab.. 
1. THE ODE 45 INTEGRATION FUNCTION 
function xdot=mateqn(t,x) 


%Define gravity 
G=386.4; 


% Define HUMMV Parameters: 

Iet=83.8; L2=46.7; 5 =05. 7 eNO a bo —294:5; 
Kel 2* 1323.0; K2=2*954.0; K3=2*1777.0; K4=2*2520.0; 
K5=1000000.0; : 
le? * 1 .2*71.5; Bit=2 12 123.05 Bee =2- 1 285.0: 
B2T=2*1.2*200.0; 


W 1=2*250.0; W2=2*250.0; W3=8200.0; 

M1=W1/G; M2=W2/G; M3=W3/G; 

11=52680; RW=5400; FW=3800; 

%Define Gun Parameters: 

L4=123.4; 5=1 1:6: LA=45.0; LS=42.6; LCG=32.5; 
L11=104.4; L12=135.0; L13=38.0; L14=17.0; L15=9.70: 
L16=110.3; 

K6=2* 1000.0; K7=2*528; KA=200000; KS=300; 
B3C=2*20.0; B3T=2* 66.0; B4=.5*250000.0; 

W4=2*200; W5=1500; W6=2300; M4=W4/G; M5=W5/G; 
M6=W6/G; 

I2=10000; I3=15000; 

Q0=12.0*pi/180; PHIA=5 1.0*pi/180; PHIS=14.0*pi/180; 


AW=3800; TB=88000; % Define Mass matrice 


%Define the mass matrix 
M=[zeros(7,7)]; 
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M(1,1)=M1;) 9 M(2,2)=M2; M(3,3)=M3+M5+M6; M(3,5)=M5*L3+M6*L3; 
M(3,6)=M5*L11+M6*L16; M(3,7)=M6*LCG*cos(Q0); 

M(4,4)=M4; = M(5,3)=M5*L3+M6*L3; M(5,5)=I11+M5*L342+M642; 
M(5,6)=M5*L3*L11+M6*L3*L16; M(5,7)=M6*L3*LCG*cos(Q0); 
M(6,3)=M5*L11+M6*L16; M(6,5)=M5*L3*L11+M6*L3*L16; 
M(6,6)=12+M5*L1142+M6*L1642+M6*L14%2; 
M(6,7)=M6*L16*LCG*cos(QQO); 

M(7,3)=M6*LCG*cos(Q0); M(7,5)=M6*L3*LCG*cos(Q0); 
M(7,6)=M6*L16*LCG*cos(Q0); 

M(7,7)=13+M6*LCG%2; 


%Define Stiffness Matrix 

K=[zeros(7,7)]; 

K(1,1)=K1+K2; K(1,3)=-K2; K(1,5)=K2*L1; 

K(2,2)=K3+K4; K(2,3)=-K4; K(2,5)=-K4*L2; 

K(3,1)=-K2;  K(3,2)=-K4; K(3,3)=K2+K4+K7; K(3,4)=-K7;  K(3,5)=- 
K2*L1+K4*L2+K7*L3; K(3,6)=K7*L12; 

K(4,3)=-K7;  K(4,4)=K6+K7; K(4,5)=-K7*L3;  K(4,6)=-K7*L12; 
K(5,1)=K2*L1;  K(5,2)=-K4*L2; K(5,3)=-K2*L1+K4*L2+K7*L3; K(5,4)=- 
Kes: 

K(9,5)=K2*01°%2+K4*12°2+K 7*L3"%2: RO .6)=khe esi 
K(6,3)=K7*L12; K(6,4)=-K7*L12; K(6,5)=K7*L3*L12; 

K (6,6)=K7*L1242+KA*LA*L13*cos(QO); 

K (6,7)=-KA*L13*LA*cos(QO); 

K(7,6)=-KA*L13*LA*cos(Q0); 

K(7,7)=+KA*L13*LA*cos(Q0O); 


%Define and velocities: 
VHF=x(10) - (L1*x(12)); VHB=x(10) + (L2*x(12)); 
VGB=x(10) + (L3*x(12)) + (L12*x(12)); 


%Define Road Inputs: 
VMPH=5.0; VEL=VMPH*17.6; PD=72/VEL; 
U1=3*(sin(2*pi*t/PD)); 


if t < L7/VEL, U2=0;else, T2=t-(L7/VEL); 
U2=3*(sin(2*pi*T2/PD)); end 


if t < L8/VEL, U3=0;else, T3=t-(L8/VEL); 
U3=3*(sin(2*pi*T3/PD)); end 


o> 


% Determine the shock absorber damping constants based on 
compression or _ tension 


% HUMMV Front shock 
if (VHF - x(8)) > 0.0, B1=BI1T; else, BI=BiC; end 


% HUMMV Rear shock 
if (VHB - x(9)) > 0.0, B2=B2T; else, B2=B2C; end 


% GUN shock 
if (VGB - x(11)) > 0.0, B3=B3T; else, B3=B3C; end 


%Define Damping Matrix 

B=[zeros(7,7)]; 

B(1,1)=B1; B(1,3)=-B1; B(1,5)=B1*L1; 

B(2,2)=B2; B(2,3)=-B2; B(2,5)=-B2*L2; 

B(3,1)=-B1; B(3,2)=-B2; B(3,3)=B1+B2+B3; B(3,4)=-B3; B(3,5)=- 
B1*L1+B2*L2+B3*L3; B(3,6)=B3*L12; 

B(4,3)=-B3; B(4,4)=B3; B(4,5)=-B3*L3; B(4,6)=-B3*L12; 
B(5,1)=B1*L1;  B(5,2)=-B2*L2;  B(5,3)=-B1*L1+B2*L2+B3*L3; 
bie )=-B3"*L3: B(5,5)=B1*L1°2+B2*L2*2+B3*L3%2; 
B(5,6)=B3*L3*L12; 

B(6,3)=B3*L12; B(6,4)=-B3*L12; B(6,5)=B3*L3*L12; 
B(6,6)=B3*L12%2+B4; 

B(6,7)=-B4; B(7,6)=-B4; —_—_B(7,7)=B4; 


% Define Forcing Vector 
F=[zeros(7,1)]; 
F(1,1)=K1*U1;  F(2,1)=K3*U2; F(4,1)=K6*U3; 


% CHECK IF FRONT HUMMV WHEEL LOSES CONTACT 
if K1*(U1 - x(1)) < -FW,  F(2,1)=-FW; = K(1,1)=K2; end 


% CHECK IF REAR HUMMV WHEEL LOSES CONTACT 
if K3*(U2 - x(2)) < -RW,  F(2,1)=-RW; -K(2,2)=K4; end 


% CHECK IF HOWITZER WHEEL LOSES CONTACT 
if K6*(U3 - x(4)) < -AW,  F(4,1)=-AW; -K(4,4)=K7; end 


%transform into system of Ist order ode's 
id=[1 zeros(1,6);0 1 zeros(1,5);0 O 1 zeros(1,4);zeros(1,3) 1 zeros(1,3)]; 
id=[id;zeros(1,4) 1 O OQ;zeros(1,5) 1 Q;zeros(1,6) 1]; 
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a=[zeros(7,7) 1d; -inv(M)*K_ -inv(M)*B]; 
fi=inv(M)* F; 
f=[[zeros(7,1)]; fi]; 


% integrate 
xdot=a*x + f; 


2. THE DRIVER CODE FOR ODE 45 


% Set the Initial Condition 
x0 1=[zeros(14,1)]; 


% Integrate in one second increments 
t01=0.00; tfi=0.999; %x0i=xn(length(xn),:); 
[ti,x 1]=ode45('mateqn',t01,tfl,x01,le-3); 
t02=1.0001; tf2=1.999; x02=xi(length(x1),:); 
[t2,x2]=ode45(‘'mateqn',t02,tf2,x02,le-3); 
t03=2.0001; tf3=2.999; x03=x2(length(x2),:); 
[t3,x3]=ode45(‘mategn',t03,tf3,x03,le-3); 
t04=3.0001; tf4=3.999; x04=x3(length(x3),:); 
[t4,x4]=ode45(‘'mateqn',t04,tf4,x04,le-3); 
tO5=4.0001; tf5=4.999; x05=x4(length(x4),:); 
[t5,x5]=ode45('mategn',t05,tf5,x05,le-3); 
t06=5.0001; tf6=5.999; x06=x5(length(xS5),:); 
[t6,x6]=ode45('mateqn',t06,tf6,x06,1e-3); 
%t07=6.0; tf7=6.999; x07=x6(length(x6),:); 
% (t7,x7J=ode45('mateqn',t07,tf7,x07, 1le-3); 
7otO8=7.0; tf8=7.999; x08=x7(length(x7),:); 
% [t8,x8]=ode45('mateqn’',t08,tf8,x08,ie-3); 


% Combine the solution state matrix 
xn=[x 1;x23x33x4;x5;x6];%x7;x8]; 
tn=[t1;t2;t3;t4;t5;t6];%t7;t8]; 

save ghl xn tn 
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APPENDIX D 
POTHOLE COURSE COMPUTER CODE 


This Appendix contains the computer code that was developed 
to model the terrain inputs for a road with a trench type pothole. This 
program calculates and plots a trace of a pothole profile and the axle 
path followed by circular wheel rolling through this type pothole. This 
program works for any pothole that can be model with a flat bottom, 
and sides that can be approximated as straight lines. The portion of 
the code shown below only works for potholes that have an entrance 
height that 1s equal to the exit height. Additional code to determine 
the terrain inputs for a any type of pothole is on file at the Naval 
Postgraduate School. Computer code was also developed to model 
bump inputs, such a log in a road and maintained on file at the Naval 
Postgraduate School. 


COMPUTER CODE: 
function [x,y,ya,u]=pothole(xent,depth,thetal ,xbot,theta2,farht,xtrail); 


% The call statement for this program/function is as %follows: 
% [x,y,ya,u]=pothole(xent,depth,thetal ,xbot, theta2,farht, xtrail) 


% Input parameters used in the call are defined as follows: 
% depth - the vertical distance from the bottom to the level 
% entrance surface. 


% thetal - the angle of the decline on the entrance side, measured 
% from an extension of the level entrance surface to the decline. 
% xbot - the horizontal distance across the bottom. 

% theta2 - the angle decline on the exit side, measured like thetal. 
% farht - the vertical distance from the bottom to the level exit 

% surface. 

% xent - the horizontal length of the level entrance region. 


% xtrail - the horizontal length of the level exit region. 


% Output parameters 


% xX - a column vector containing a discrete amount of evenly spaced 
%o values that were used to calculate the profile and axle path. 
7 y - a column vector that contains the ‘y' values for the pothole 
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vy profile. The bottom of the pothole is defined as "zero" y value. 


To ya - a column vector containing the axle path 'y' values. 
To u - a column vector containing the height of the bottom of the 
Yo wheel. 'u' is the input for the simulation model 


% Other parameters assigned in the program: 
Jo deltax - the ‘x’ intraval spacing. 


% xent  - the horizontal length of the level entrance region. 

To xtrail - the horizontal length of the level exit region. 

% rad - the radius of the wheel that is rolling through the pothole. 
To Varibles that are calculated based on the above assigned 

% parameters: 

To 

7o rhtetal & rtheta2 - the thetal & theta2 in radians. 

Yo xnear & xfar - the x distance differance from the top to the 
%o bottom for the near & far sides. 

7o xcritl & xcrit2 - the x value where the wheel can no longer 
To follow a circular path as it goes around the near & far 

% corners(near into hole, far out of hole). 

To xcrl - the x value where the wheel hits the far side or corner 
% xcr2 - the x value from where the wheel hits the the exit region 
To from the near side or near corner. 

To xchl - the x value at the bottom corner on the entrance side. 
%o xch2 - the x value at the bottom corner on the exit side. 

%o xch3 - the x value at the top corner on the exit side. 

7o xend - the ending x value of the pothole profile. 

Zo ycornl & ycorn2 - the y values when the axle is following a 
7o circular path around the near & far corners respectively 


Jo *¥*KKEEKaSSIOnMment of pothole & wheel parameterst*********"**4aus 
deltax=.01; rad=17.0; 


q 3K 3K 2K 2K 2K 3K ok ok ok calculating other varibles used 3K AK 2 OC oie I OC ie aK OK ie 2 OI OC 2k 2K 2 i 2K 2 2k 2K 2K 2k 


rtheta1=thetal *pi/180;, rtheta2=theta2*pi/180; 
xnear=depth/tan(rthetal);, xfar=farht/tan(rtheta2); 
xchl=xnear+xent;, xch2=xch1+xbot; 

xch3=xch2+xfar;, xend=xch3+xtrail; 
xcritl=xenttrad*sin(rthetal);, xcrit2=xch3-rad*sin(rtheta2); 
xcrit3=xent+trad*sin(rtheta2);, xcrit4=xch3-rad*sin(rtheta2); 


% converting/rounding x-switch values to integer increments of 
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% deltax; keeps the output x, y,&ya vectors the same size 


num=xent/deltax;, num=round(num);, xent=num* deltax; 
num 1=xch1/deltax;,num1l=round(num1);, xchl=num |! *deltax; 
num2=xch2/deltax;,num2=round(num2);, xch2=num2*deltax; 
num3=xch3/deltax;,num3=round(num3);, xch3=num3*deltax; 
num4=xend/deltax;,num4=round(num4);, xend=num4*deltax; 
xdl=xcritl-xent;,nl=xd1/deltax;,nl=round(n1);, xcritl=xent + 
n1l*deltax; 

xd2=xch3-xcrit2;,n2=xd2/deltax;,n2=round(n2);, xcrit2=xch3 - 
n2*deltax; 


(a. ~*~ creating x vectors for pothole regions *************exx* 


x 1=0:deltax:xent;, sizl=size(x 1); 
x2=xent+deltax:deltax:xchl;, 
x3=xch1+deltax:deltax:xch2;, §1z3=size(x3); 
x4=xch2+deltax:deltax:xch3;, 
x5=xch3+deltax:deltax:xend;, §1zZ5=size(x5); 


x=[xl x2 x3 x4 x5]; 


% xe Calculating y values for the pothole profile **********xxxx*x 
y l=depth*ones(1:siz1(2)); 

y2= -tan(rthetal)*x2 +tan(rthetal)*xent + depth; 

y 3=zeros(1:s1z3(2)); 

y4= tan(rtheta2)*x4 - tan(rtheta2)*xch2_ ; 

y5S= _ farht*ones(1:siz5(2)); 

y=lyl y2 y3 y4 y5]; 


% ****** calculating the axle trace in the entrance region******** 
xal=0:deltax:xent;,sizal=size(xal); 
yal=(depth+rad)*ones(1:sizal(2)); 

im xcrit2 <= xcritl, 
%***checking if the axle motion is only circular***** 
if xcrit2 > xent, 
%finding the x intercept of the near & far side circular motion 
for xint=xcrit2:deltax:xch3, 
yl=depth + sqrt((rad.*2) - ((xint-xent).*2)); 
y2=farht + sqrt((rad.*2) - ((xint-xch3).*2)); 
if yl<=y2; break,end 
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end 


else 
for xint=xent+deltax:deltax:xch3, 
yl=depth + sqrt((rad.42) - ((xint-xent).*2)); 
y2=farht + sqrt((rad.“2) - ((xint-xch3).%2)); 
if yl<=y2; break,end 
end 
end 


% calculating the near and far side y values for the axle path 
xa2=xent+deltax:deltax:xint; 
ya2=depth + sqrt((rad.*2) -((xa2-xent).%2)); 
xa3=xint+deltax:deltax:xch3; 
ya3=farht + sqrt((rad.“2) -((xa3-xch3).*2)); 
yac=[ya2 ya3] ; 
xa4=xch3+deltax:deltax:xend;, § siza4=size(xa4); 
ya4=(farht+rad)*ones(1:siza4(2)); 
ya=[yal yac ya4]; 
%checking & calculating axle path when its not just circular 
else 
7o****intializes a counter & the vector yac**** 
n=0;, nnum=(xch3-xent)/deltax;, yac=zeros(1,nnum); 
7o*****checks & determination of axle path**** 
for xa=xent+deltax:deltax:xch3, 
if xa <= xcritl, 
ycornl=depth + sqrt((rad.*2) - ((xa-xent).*2)); 
yfar=tan(rtheta2)*xa-tan(rtheta2)*xch2+rad/cos(rtheta2); 
if ycornl > rad & ycornl > yfar, yacl=ycorn]; 
elseif yfar > rad, yacl=yfar; 
else yacl=rad; end 
elseif xa >= xcrit2, 
ycorn2=farht + sqrt((rad..2) - ((xa-xch3).*2)); 
ynear=-tan(rthetal)*xa+tan(rtheta1)*xent... 
+depth+rad/cos(rtheta1); 
if ycom2 > rad & ycorn2 > ynear, yacl=ycorn?; 
elseif ynear > rad, yacl=ynear; 
else yacl=rad; end 
else 
ynear=-tan(rtheta1)*xa+tan(rthetal)*xent... 
+depth+rad/cos(rtheta1); 
yfar=tan(rtheta2)*xa-tan(rtheta2)*xch2+rad/cos(rtheta2); 
if ynear > yfar & ynear > rad, yacl=ynear; 
elseif yfar > rad, yacl=yfar; 


el 


else yacl=rad; end 
end 
(lor) Lae yac(1,n)=yac 1; 
end 
Oy HF Ae A He 8 He Ae Ae ae He 2 Ae Ae axle y values in the exit region 3K OK 2K OK 2K OK ok OK 3K OK OK OK ok Ok 
xa4=xch3+deltax:deltax:xend;, $1za4=size(xa4); 
ya4=(farht+rad)*ones(1:siza4(2)); 
ya=[yal yac ya4]; 
end 
Jo TEXKEKKKEE DT OT OF POTHOLE PROFILE & AXLE PATH ********#*4#% 


ax=[0,xend,-50,yplot]; 

y=y-depth; ya=ya-depth;u=ya-17; 
X=X';y=y 5ya=ya';u=u; 

matrix=[x y ya ul]; 

save road2dat.m matrix 

ax1s(ax); 

plot(x,y); hold on ; plot(x,ya) ,pause,hold off 


APPENDIX E 
NATURAL FREQUENCIES AND MODE SHAPES 


Table E-1 below contains the natural frequencies and mode 
Shapes of the HUMMV/Howitzer system. They were calculated using 
the procedures decsribed in Chapter IV. In Table E-1 the system's 
natural frequencies are listed in the top row. The _ positional 
coordinate values of the mode shape which corresponds to a 
particular natural frequency are listed in the column below the 
natural frequency. The natural frequencies shown in Table E-1 have 
units of hertz. The positional coordinates have units of inches for Z), 


Z,,Z3, and Z,, and radians for the rotational cordinates 6;, 02 and 63. 


TABLE E-1 NATURAL FREQUENCIES AND MODE SHAPES 


Natural Frequencies 





coe 1.286 1,870 3,328 12.26 13.61 18.51 2one0 
Zi -0001  .0783 -.2392 -.0212 1.217 + .0130 ~—.0342 
Zo 0130 .1293 1425. 0076 = 0364S 1.222. -. 1221 
Z 0140 .2057 -.0408 0018 -.0221 .0138 0277 
Zh, 1064 0168  .0202 -1.385 -0.221 -.0097 0341 
Q, 0002 .0003 .0059 .0001 .0011 0005 .0025 
Q, -0025 -.0013 -.0026 .0001 ~-.0003 -.0001 -.0036 
03 -0025  -.0013 -.0027 .0003 -.0006 -.0011  .0063 


APPENDIX F 


A-FRAME CALIBRATION DATA 


The following tables contain the load versus strain data that 
was recorded during the A-Frame calibration procedures that are 
described in Chapter V. The locations of the strain gages is shown on 
Figures 5-1 and 5-6. Tables F-1 and F-2 contain the data from the 
vertical load tests and Tables F-3 and F-4 contain the data from the 
lateral load tests. The method used to calculate stress yalues from 


the strain gage readings is described in Chapter V. 


TABLE F-1 VERTICAL TEST #1 A-FRAME CALIBRATION DATA 
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TABLE F-1 VERTICAL TEST #1 A-FRAME CALIBRATION DATA (CONT) 
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TABLE F-3 LATERAL TEST #1 A-FRAME CALIBRATION DATA (CONT) 





a 
= 
4 
o 
~ 
lJ 
g 
< 
© 


UW 
= oe 
~ 
4 S 
~) aay 
~) 
oD < = 
te 2 
= 
(¢ 9) - isp) 
= a ‘t 
< id] © 
< < 
© os) ©) 


FORCE (LBS 


120 


TABLE F-4 LATERAL 


-649 


-379 | -530_ |-643 | -539 | -446 | -364 | -276 


-745 


A-FRAME LATERAL LOAD TEST 2 
-276 


975 


117 





FORCE (LBS 
AVG STRAIN 1&2 


4.46 


TEST #2 _A-FRAME CALIBRATION DATA 


0.96 


~ 
od es ©] ©] — 
=i = =| oO 
co @ 
= ~ 
QJ) © 
is 
en LB 
1 o> 
O! Oo] © || Se 
©} ©} i 


4.66 


0.79 
_ 
Eee 





ell I ad | 
Oy} wor wo ~~) 
EOL 5 WH is 
wz wo] © Ol oO] © 
OS) = So ied Ld (if 
rN] ©) © © Ol wo Ol olo 
wT 
© , 
ro oO re) 
on ] 09] 9 
foo] aa’ we © 
mS] @ wo] 2 = 
Ql = NN — 
Oi! © 
| = 
w 
= 1“ or) = ¢ = 
7) 7) 7) 7) 
ax |S a|<} |S ax} |S o| 3} ¥ 
a col x i 4 =| 
wo wo 7p) w 
oie) | fe vel dal | ta del dal | fee 
” Vo} ep) ola Ow 
wae! {2 zjH El 12 zjA1~} 12 Z| O15 
LU < —=|W) w < =| Ww < az| Win 
elo} (2 TEI) 1° Z\E}o} |o IE 
OZ |Z Fil 2Z| [Zolli B72) [= | |Z 
all |Z 2Z}S) |<jujw! |sjo} i< z|9 
<jz c ol </ 2 ra oi <2 rd | <| oa 
><} Wy (= >] ><} WW fe ellie ~ >} ><] Ww 
<iao} in <i</ao] jo <j<cio| |w <| <|/o 





fee 


TABLE F-4 LATERAL TEST #2 A-FRAME CALIBRATION DATA (CONT 


N 12) 
ald 


S 
w 
~~ 


SACTTPH rg ITA ALL 


RANSVERSE 





GAGE 4 saa 


Y) 
02) 
Lu 
a 
Ke 
9) 
a 
< 
a 
2 
a 
© 


FORCE (LBS 
GAGE 3 (CLAMP 
CLAMP STRESS (KSI 
GAGE 11 

GAGE 13 (AXIAL 


122 


APPENDIX G 


IOWITZER ORIENTATION PICTURES AND DIAGRAMS 


This Appendix contains pictures and diagrams that identify the 


key components of the M119 Howitzer that were important for this 


Study. 





Figure G-! Towing in Firing Position 
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Figure G-3 A-Frame/T-Bar Support 








Figure G-4 T-Bar (View from top Down) 
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Figure G-5 
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